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ABSTRACT 


Low-frequency variability is examined with a focus on the 30-60 day oscillations in 
General Circulation Model (GCM) integrations made with the GLAS GCM, the UCLA 
GCM, and the GLA GCM. A comparison with global ECMWF analyses of 1984-1986 
was carried out. One purpose of this study is to illuminate some key issues critical to 
present theories about intraseasonal oscillation through intercomparison of the three GCMs 
and observations, and thus to provide insight for improving the critical parametenzations 
needed for better simulation of low-frequency variability. All three GCMs are grid point 
models with a 4 x 5 degree latitude/longitude resolution and with 9 layers in the vertical. 
The models were integrated for at least two years in the seasonal cycle mode. Each of the 
three GCMs employs a different scheme for computing finite differences. 

The GLAS GCM simulates a realistic eastward propagation of the 30-60 day 
oscillation in the tropical upper-level divergent flow. The oscillation becomes quasi- 
stationary over the Indonesian region and accelerates over the central Pacific, as observed. 
In the GLA GCM, the oscillation propagates into the higher latitudes of both hemispheres 
as the waves leave the convective region. The eastward propagating oscillation is not 
obvious in the UCLA GCM. The wavenumber-frequency spectra of 200mb velocity 
potential reveal that all the GCMs have a significantly weaker signal of eastward 
propagation in the 30-60 day range than do the analyses. The spectrum for the GLAS 
GCM is dominated by 20-60 day periods, while the GLA GCM has a spectral peak around 
the 20-30 day periods. There is a weak peak near 15 days for eastward propagation in the 
UCLA GCM. The relationship between dominant phase speeds of the oscillation and the 
vertical profile of heating in the GCMs is in general agreement with current theory. 

The composite patterns of the observations indicate that a Kelvin wave-type 
structure is dominant in the tropics near the center of the oscillation. The simulated winds 
are fairly realistic, although the meridional component is too strong, especially in the GLA 
GCM. The vertical structures of the zonal component of wind and moisture suggest that 
mobile wave-CISK (Lau and Peng, 1987) is an important mechanism in maintaining the 
intraseasonal oscillation in these GCMs. The vertical distribution of the moisture field 
further suggests that evaporation-wind feedback (Neelin, et al., 1987) may play a role in 
maintaining the eastward propagating waves in the tropics. The differences in the 
structures of the oscillation in the GLAS GCM and GLA GCM appear to be a consequence 
of the different numerical schemes used. The GCMs have preferred zones for diabatic 
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heating, with turn-on heating when the rising branch of the intraseasonal oscillation passes 
over these convective regions. 

All three GCMs fail to capture the detailed evolution of the different stages of the 
development and decay of the oscillation. Results suggest that an improvement of the 
boundary layer moisture processes may be crucial to a better simulation of the oscillation. 
In addition, both the deficiencies in basic state circulations of the models and deficiencies in 
the 30-60 day oscillation appear to produce the unrealistic tropical-extratropical interactions 
in the GCMs. 
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I. INTRODUCTION 


in recen, years, .here have been a large number of GCM simulation experiments of 
tropical intraseasonal oscillations. In general, these experiments fall into two categones. 
First, are experiments based on realistic boundary conditions such as observed sea surface 
temperatures, snow-ice distributions and topographic forcing (e.g. Lau and Lau, 1986, 
Hayashi and Colder, 1986; Geisler and Pitcher, 1988; Von Storch et al„ 1987, Tokioka e. 
ai„ 1988 and others). Many of the experiments in this category were not initially designed 
for studies of intraseasonal oscillations per sc, but were motivated by more general climate 
concerns. Second, are simplif.ed GCM experiments which were designed to study the 
basic mechanisms of intraseasonal oscillations using idealized boundary conditions 
aqua-planet experiments (Hayashi and Sumi, 1986, Swinbank e, al„ 1988 ) and swamp 
ocean experiments (Lauetal., 1988). Both categories of experiments have provided rich 
sources of information and testing grounds for validating theories of intraseasonal 
oscillations (e.g. Lau and Peng, 1987, Wang, 1988, Chang and Lim, 1988, Emanuel, 


1987, Neelin et al., 1987, Hendon, 1988 and many others ). 

The results of die above diverse GCM experiments are generally consistent with the 
simple idea that the structure and eastward propagation of inoaseasonal oscillations are due 
to the interaction between low-level moisture convergence and mid-tropospheric latent heat 
release from deep convection, i.e., the so-called mobile wave-CISK mechanism (Lindzen, 
1974; Lau and Peng, 1987). However, a number of discrepancies still exist between 
model results and observations. One of the most serious difficulties is that almost all the 
GCM results produce oscillations with associated periods which are about 20 50V 
than the observed results (20-30 days instead of 30-60 days). Lau and Peng (1987), Sui 
and Lau (1989) anil Takahashi (1987) suggest that the phase speed of propagation is 
sensitive to the level of maximum heating and that a lower level of heating will give rise to a 
slower phase speed of the oscillations. More recently Lau and Shen (1988) suggest that 
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some degree of coupling with [he ocean may be required to produce the observed phase 
speed of these oscillations. On the other hand, Hayashi and Colder (1986) indicate that a 
distinct low-frequency signal at 40 days can be obtained by simply increasing the spatial 
resolution of the GCM. Overall, it appears that all COM simulations indicate that some 
form of mobile wave-CISK mechanism is at work. Yet this mechanism is unable to 
explain the observed continual eastward propagation around the globe of the upper level 
wind field which occurs in spite of the lack of latent heal forcing over the eastern Pacific 
sector. Another problem is that the most -realistic" simulation of the intraseasonal 
oscillations appears to be associated with a tum-on latent heating over the maritime 
continent and tropical western Pacific rather than a propagating heat source as observed 
over these regions. The aqua-planet experiments simulate very well the combined eastward 
propagation of convection and circulation but fail to reproduce the stationary component of 

the oscillations. In reality, both propagating and stationaty components are present in these 
oscillations (Lau and Chan, 1985). 

Obviously, much work is still needed in GCM modeling, theory and observations 
in order to fully understand the dynamics of intraseasonal oscillations. In view of the 
diverse model structure and physical parameterizations in different GCMs, there is a need 
to compare the results of the c l imato l ogy of intraseasonal oscillations in various GCMs as 
they relate to the model characteristics. This will provide us with clues for the model 
differences and hopefully, new insights for improving the critical parameterizations needed 
for better simulation of these oscillations. Archives of multiyear GCM climate integrations 
that now reside in various research laboratories can provide a rich repositoiy of information 
from which these intercomparison studies can now be conducted without further straining 
resources needed for extended integrations. In this paper, a modest step is taken in this 
direction. Obviously, the results of this work are not going to resolve in any definitive way 
the aforementioned problems, but they will illuminate some key issues critical to present 
theories of intraseasonal oscillation and in doing so, demonstrate a methodology by which 
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model intercomparison can be carried out in the future. For this paper, we have had to 
work under the constraints that only a limited set of common parameters had been saved 
and that due to institutional changes and other factors, the integrations cannot be repeated 
exactly or extended to pin down the causes of specific differences. However, based on the 
present results, we shall provide some guidelines for further study of GCM intraseasonal 
oscillations. As our personal preference, we shall use the general terminology 
"intraseasonal oscillation" to refer to the 30-60 day or 40-50 day oscillation, or the Madden 
and Julian Oscillation, names which are also commonly used in the literature. They will be 
used interchangeably in the paper. The seasonal mean patterns, latitude-pressure sections 
and longitude-pressure sections for selected parameters are shown in the Appendix. 
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II. DESCRIPTION OF THE GCMS AND OBSERVATIONAL DATA 


This study examines low frequency variability in General Circulation Model (GCM) 
integrations with the GLAS GCM, the UCLA GCM, and the GLA GCM. This is compared 
with global ECMWF analyses of 1984-1986. All three GCMs are grid point models with 
a 4 x 5 degree latitude/longitude resolution and with 9 layers in the vertical, and were 
integrated for at least two years (three in the case of the UCLA GCM) in the seasonal cycle 
mode. That is a smooth seasonal cycle in insolation including diurnal variation and all the 
other boundary conditions were introduced into the model integrations. The GLAS and 
GLA models use identical formulations of all physical processes, including radiation, 
boundary layer processes and cumulus convection. The parameterization of physical 
processes in the UCLA GCM are completely different and are substantially more 
sophisticated. Each of the three GCMs employs a different scheme for computing finite 
differences. They all use the vertical differencing of Arakawa and Suarez (1983). 

The GLAS GCM is described in some detail by Randall (1982). This model 
employs a staggered grid (the so-called B grid), and the second order finite-difference 
scheme for horizontal momentum advection conserves momentum and kinetic energy. The 
cumulus parameterization scheme of Arakawa (1969, 1972), developed originally for the 
three-level UCLA GCM, was modified for use in a nine-level model by Somerville et al. 
(1974). This was accomplished by combining the lowest six layers pairwise into three 
layers for the purpose of cumulus convection, and restricting this process from extending 
into the uppermost three layers. This constraint (also present in the GLA GCM, to be 
described below) strongly influences the simulation of intraseasonal oscillations in the 
tropics. The boundary layer parameterization follows Deardorff (1972) with the 
assumption that the depth of the boundary layer is 500 meters everywhere. Vertical 
diffusion of momentum, sensible heat, and moisture above the boundary layer are 
negligible. Long-wave radiation (Krishnamurthy, 1982) was called once every five hours, 


4 



but the heating rates were applied every 30 minutes, while short-wave radiation (Davies, 
1982) was called every 30 minutes. The two-year integration studied in this paper was 

described in detail by Straus and Shukla (1988a, 1988b). 

The GLA GCM is identical to the GLAS GCM in terms of all physical 
parameterizations, and so differs from it only in terms of the numerical methods used. The 
GLA GCM utilizes an unstaggered grid (the so-called A grid) with a quadratic (energy) 
conserving scheme of fourth-order accuracy (Kalnay « <rU983). Sixteenth-order Shapiro 
filters are used in the horizontal in both the GLA and GLAS GCMs, although they are used 
in different ways. With the staggered grid (GLAS GCM) the filter is not necessary to 
integrate the model, but the filter is necessaty to prevent a divergence of the solution when 
the unstaggered grid (GLA GCM) is used (Randall, personal communication). 

The UCLA GCM (described in detail in Randall, et al., 1985) uses a more recent 
parameterization of cumulus convection due to Arakawa and Schubert (1974) and Lord et 
al. (1982). One distinctive feature of this GCM is that the planetary boundary layer top is a 
prognostic variable and acts as a coordinate surface. Mass sources and sinks for the 
planetary boundary layer are large-scale divergence, turbulent entrainment and cumulus 
mass flux (Randall, 1979, 1984). The surface flux determination is that due to Deardorff 
(1972), and the radiatton that of Schlesinger (1976). The horizontal differencing scheme 
uses a staggered grid, the so-called C grid. In this scheme, nondivergent flow conserves 
kinetic energy, while potential enstrophy is conserved by the full flow. In addition, 

advection of vorricity and temperature are of fourth-order accuracy. 

While the global analyses prepared by the European Centre for Medium-Range 
Weather Forecasts (ECMWF) since 1980 were available, we chose to use only the more 
recent analysis of 1984-1986 because prior to 1984 the interannual variability of the 
divergent wtnd component was no. reliable (Trenberth and Olson, 1988). In this study we 
worked primarily with five-day average fields for both the GCMs and the observations 
(analyses). In each case the time mean (over the entire length of the dataset) and the annual 
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and semiannual harmonics were removed. (These harmonics were calculated by least- 
squares fitting for each year separately, followed by averaging over the individual years.) 
For many of the diagnostics presented, the time series of five-day means was further 
filtered by calculating the raw Fourier components, retaining only the periods between 30 
and 60 days, and performing the inverse transform. 
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III. GENERAL FEATURES OF INTRASEASONAL VARIABILITY 


a. 200 mb velocity potential 

Following prev ious work, we use the 200 mb velocity potential as the primary field 
with which to detect the presence of the tropical intraseasonal oscillations. Figures 1 and 2 
show the geographical distribution of the total variance of this field (after the seasonal cycle 
is removed) and the low-frequency variance (filtered to retain only periods of 30-60 days), 
respectively for the three GCMs and the ECMWF analyses. The analyses [Figures 1(d) and 
2(d)] show that the total variance maxima occur in the Indian Ocean and western Pacific 
region and that the strongest degree of 30-to-60 day variability occurs over the Indian 
Ocean, with large variability also seen over the western Pacific, and extending across the 
equatorial Atlantic from South America to western Africa. This pattern is similar to the 
structure found in the observational studies of Lau and Chan (1985, 1986) using outgoing 
long-wave radiation iOLR). The GLAS GCM results [Figures 1(a) and 2(a)] indicate 
enhanced variability ever the western Pacific (although the maximum is shifted about 30 
degrees too far east) and over South America. Noticeable is the absence of a distinct 
maximum in variability over the Indian Ocean, the region where the observed signal is the 
largest. The magnitude of the simulated variability is quite weak compared to observations. 
This is even more so in the case of the UCLA results [Figures 1(b) and 2(b)]. The GLA 
GCM [Figures 1(c) and 2(c)] is most realistic in terms of the overall variance pattern, but 
the simulated variance is still not as strong as is observed, especially over the Indian Ocean. 
The simulated maxima in the western Pacific region in the GLAS and GLA models are 
closely tied to the locations of the maxima in mean diabatic heating in these GCMs. 

The variance distribution in the frequency domain is examined by spectral analysis. 
In Figures 3-6, each panel shows the spectral density of 200 mb velocity potential at the 
location represented by longitude and latitude (dotted) lines. The analyses [Figure 3] show 
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that the spectral power is pronounced in the period range of 20-70 days, with double peaks 
around 30 day and 60 day periods, in the extensive region over the Indian Ocean and 
western Pacific. These spectral peaks are substantially weaker in the central Pacific and 
high latitudes except over East Asia. The GLAS GCM [Figure 4] shows a realistic spectral 
distribution but much weaker 30-60 day oscillation (note different scales of y-axis). The 
spectral peaks are shifted toward higher frequencies. In the UCLA GCM [Figure 5], the 
spectra have no concentrated power in the 30-60 day period range. In the GLA GCM 
[Figure 6], the 30-60 day oscillation is pronounced only in the limited region over the 
western Pacific. 

The presence of eastward propagation is described by a Hovmoller diagram of the 
five day mean 200 mb velocity potential in the tropics (20°S-20°N) in Figures 7 - 10 for 
the GCMs and the ECMWF analyses. The analyses [Figure 7] shows more coherent 
eastward propagation during the northern winter. The oscillation is enhanced in the Indian 
Ocean and western Pacific region and in the Atlantic. The GCMs [Figures 8-10] show 
eastward propagation, but the amplitude is substantially weaker. The GLAS GCM [Figure 
8] properly simulates the enhanced amplitude and more coherent propagation during the 
northern winter. As indicated earlier, the UCLA GCM [Figure 9] shows very weak 
amplitude of eastward propagation. In the GLA GCM [Figure 10], the oscillation has too 
strong a standing component in the Indian Ocean and western Pacific and it is not 
weakened during the northern summer. 

The eastward propagation feature is also revealed by one-point correlation maps of 
the five day mean 200 mb velocity potential. The correlations are computed from the 
filtered data with respect to a base point at the Equator and 120°E, as shown for the GCMs 
and the ECMWF analyses in Figures 11-14. Lags of - 10 to + 10 days are considered. 
The propagation is quite evident in Figure 1 1 for the observations, as the planetaiy scale 
(zonal wavenumber one) structure gradually moves eastward, appears to become more 
stationary over the Indonesian region, and then speeds up over the central Pacific. This is 
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consistent with the observed propagation of OLR and wind patterns for the 30 - 60 day 
oscillation (e.g. Lau and Chan, 1985; Knutson et al., 1986). Both theoretical analysis (Lau 
and Shen, 1988) and observational study (Gutzler and Madden, 1989) have suggested that 
the reduced phase speed of the disturbance may be related to moist Kelvin wave modes 
associated with enhanced convection over the warm SSTs in the western Pacific, and the 
increased phase speed further east to dry circulation modes due to the lack of convection 
over the cooler waters. 

Figure 12 indicates that the GLAS GCM successfully simulates the eastward 
moving 30 - 60 day oscillation in the upper level divergent flow, including the acceleration 
of the wave after it passes the convectively active region. The one-point correlation maps of 
the UCLA GCM (Figure 13) are much less realistic. Only local stationary signals are 
apparent. These signals disappear at ±10 days. For the GLA GCM (Figure 14), the 
waves are nearly stationary at the base point in the tropics, but propagate eastward 
dispersively into the higher latitudes of both hemispheres as the wave leaves the Indonesian 
region. The zonal wavenumber one character, so evident in the correlation maps of the 
ECMWF analyses, is present in the results from the GLAS and GLA GCMs, but is much 
less obvious in the UCLA correlations. 

The dominant space and time scales of the eastward traveling waves can be further 
elucidated by the use of wavenumber - frequency spectra for the 200 mb velocity potential. 
Figures 15-17 show the wavenumber-frequency spectra at the Equator, 30°N and 30°S, 
respectively. The spectra were calculated using a Hamming window with a truncation 
point equal to half the total number of data points. In these figures, the ECMWF spectrum 
is dominated by zonal wavenumber one contribution and eastward propagation with 
periods of 30 to 60 days, with peaks at 30 and 60 days. The spectrum for the GLAS GCM 
is dominated by similar time-scales (20 to 60 days) for eastward propagation, but the 
overall strength of the signal is significantly weaker than observed, consistent with Figure 
2. There is no indication of a distinct oscillation in the 30 to 60 day range in the UCLA 
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GCM spectrum, although an eastward propagating peak near 15 days is suggested at the 
Equator and 30°N. The GLA GCM, which has more power in the 30 to 60 day range than 
does the GLAS GCM but less than the ECMWF analyses, is dominated by eastward 
moving waves with a peak around 20 - 30 days. Thus while the GLA spectrum shows 

more power in the 30 - 60 day range than the GLAS spectrum, the overall shape of the 
latter is more realistic. 

The latitude-frequency distribution of spectral power for zonal wavenumber one in 
the ECMWF analyses [Figure 18] shows that the 30-60 day oscillation is global in extent 
and nearly symmetric about the Equator. Planetary-scale waves in the 200 mb divergence 
field with 30-50 day periods were also described by Lorenc (1984) using ECMWF 
analyses of FGGE year. However, any equatorially trapped divergence field may possess 
a much larger scale velocity potential because of the relationship between the two fields 
(Hendon, 1986). It is noted that the eastward propagating 30-60 day waves in the velocity 
potential of the GLAS GCM and the UCLA GCM are more closely confined to the Equator 
than is the case in the analyses. The latitude-frequency spectra for zonal wavenumber two 
in the ECMWF analyses [Figure 19] show the dominant eastward propagating waves 
around 45-60 day period range in the Northern Hemisphere. It is noted that the GLA GCM 
has unrealistically strong eastward propagating waves in the Southern Hemisphere, with 
periods around 45 days and 23 days at 40°S and period around 11 days at 30°S. The 
latitude-frequency spectra for the zonal wavenumbers 3-10 [Figure 20] show much weaker 
signal in the 30-60 day oscillation in the tropics with a comparable power in the eastward 
and westward components. The GCMs and analyses show somewhat stronger eastward 
travelling component off the Equator in the Southern Hemisphere. 

b. Vertical heating profile 

In a series of simple numerical experiments, Lau and Peng (1987) and Sui and Lau 
(1989) showed that the phase speed of the intraseasonal oscillation depends in a sensitive 
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way upon the vertical structure of the convective heating. In order to define a diabatic 
heating associated with the oscillation in the GCMs, we composited the total diabatic 
heating on the basis of the time series of the filtered (five-day mean) velocity potential 
averaged over the western tropical Pacific (10’S - 10°N, 120°E - 180°), and included in 
the composite only periods in which the velocity potential exceeded one standard deviation 
below the sample mean. (A minimum in the velocity potential is associated with 
divergence.) Longitude/height sections of the composite total diabatic heating are shown 
for all three GCMs in Figure 21. While the GLAS and GLA GCMs show much stronger 
heating rates in the composite than does the UCLA GCM, there is an important difference 
in vertical structure as well. The heating profile of the UCLA GCM is nearly uniform with 
height throughout muc h of the troposphere. In contrast, the other two models have their 
heating concentrated in the mid-troposphere, with a well-defined maximum near 500 mb 
and small heating rates at upper levels. 1 However, the vertical profile of composite 
heating averaged over the above-mentioned region [Figure 22] reveals that the region of 
largest heating in the GLAS GCM heating extends further downward than in the GLA 
GCM, so that while the latter model produces the stronger heating at its peak (near 450 
mb), the GLAS GCM heating is stronger at 700 mb. Furthermore, the GLAS GCM 
heating exhibits only one center while the GLA GCM shows double maxima separated by 


30° longitude. 

According to a simple theory by Sui and Lau (1989), the downward displacement 
of the GLAS GCM s heating profile in comparison to that of the GLA GCM should be 
reflected in the phase speed of the wave, with the latter (GLA) model having the faster 
phase speed. This is clearly seen in the spectra, as discussed above. The absence of an 
intraseasonal oscillation signal in the 30 - 60 day range in the UCLA GCM is consistent 
with the nearly uniform profile of vertical heating, as shown by Tokioka et al. (1988) in the 


iThe form of the vertical profile of diabatic heating in the GLAS and GLA GCMs is a dircct “" se ^ UC ° 
the manner in which the Arakawa cumulus convection scheme was implemented in these models. I 
particular, no convection was allowed to penetrate into the top three model layers. 
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context of experiments with a five-level aqua-planet model utilizing the same Arakawa- 
Schubert cumulus convection scheme as in the UCLA GCM. These authors find only a fast 
eastward propagating wave (with a period of 10 days), which is to be compared with the 
period of 15 days we find in the UCLA model. However, by setting a lower limit to the 
entrainment of environmental air in the cumulus parameterization, thereby inhibiting very 
deep convection, Tokioka et al. show that a 30-60 day oscillation then does appear in 
conjunction with diabatic heating profiles which have a distinct peak in mid-troposphere. 

The longitudinal positioning of the largest heating in the composites for each of the 
three GCMs coincides with the maximum of velocity potential variance in the western 
Pacific for that GCM seen in Figure 2. This positioning also coincides with the longitude 
of maximum mean equatorial heating in each model, as seen from maps of seasonal mean 
precipitation [Appendix A25-A27]. This strongly suggests that the diabatic heating 
associated with the intraseasonal oscillation is triggered as the sector of the wave with low- 
level convergence passes the region in which the model tends to produce heating in any 
case. Similar turn-on heating is found in a number of simulations of the 30 -60 day 
oscillation by GCMs which incorporate realistic lower boundaries (e.g. Lau and Lau, 1986; 
Geisler and Pitcher, 1988; Von Storch, et al. 1987). This should be contrasted to results 
of aqua-planet models (i.e., those with axisymmetric ocean or swamp lower boundaries), 
in which the precipitation follows the wave (Hayashi and Sumi, 1986; Swinbank et al. 
1988; Lau, et al., 1988). The observational study of Lau and Chan (1985) showed 

evidence of both propagation and stationary components of convection associated with the 
30-60 day oscillations. 
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IV. COMPOSITE ANALYSIS 


In order to examine the detailed structure of the eastward propagating oscillation as 
it appears in the Indian Ocean and western Pacific sectors of the models and in the 
analyses, we use a compositing scheme which portrays the oscillation in a frame of 
reference moving with the wave. This scheme is a modified form of the one introduced by 
Lau et al. (1988), who studied this oscillation in a model with axisymmetric boundary 
conditions and so did not confine their attention to one particular sector. The first step of 
this scheme is to average each five-day mean velocity potential (from the filtered series) 
from 10°S to 10°N and retain only the first three zonal wavenumbers, yielding a smooth 
field which varies only in longitude. A translation T~ is defined which maps the minimum 
of the smoothed velocity potential to the origin of longitude. A second translation T+, 
which maps the maximum to the origin, is also defined. The translations T" and T + are 
applied to whatever field (F) is to be composited, leading to the translated fields F~ and F^ , 
respectively. The difference (F‘ - F + ) is taken and is averaged over all available pentads, 
regardless of season. In this average, each pentad is weighted by the longitudinal rms 
deviation of the field F, thus serving to emphasize periods when the oscillation is strong. 
This average weighted difference defines the oscillation in a frame of reference attached to 
the moving wave. The only modification we have introduced to the scheme of Lau et al. is 
to exclude from the average those periods when the minimum of the velocity potential did 
not fall in the sector 60°E - 180°, so that we are looking at the structure of the oscillation 
only in the Indian Ocean - western Pacific sector. 

This series of steps has been repeated using temporal lags or leads, so that the data 
fields being composited and the longitude shifts used refer to different times. More 
specifically, only one longitudinal coordinate system is used for the various time lags of a 
given composite, and that system gives the position with respect to the arrival of the 
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composite wave at zero lag. The non-zero lag composites show the fields some days before 
or after the arrival of the wave in this wave coordinate system. Further details are given in 
Lau et al. (1988). 

a. Horizontal structure of wind and height 

The upper level (200 mb) composite tropical circulation is shown for the GCMs and 
the analyses at zero lag in Figure 23. The reader should keep in mind that the composite 
essentially follows the wave within the broad Indian Ocean - western Pacific region, so that 
longitude in Figure 23 and subsequent figures is a coordinate with respect to the center of 
the wave. The ECMWF analyses show a very broad region of outflow dominated by the 
zonal component of the wind, with easterlies to the west and westerlies to the east of the 
wave center. The structure near the Equator clearly resembles that of a Kelvin wave, with a 
possibly weak Rossby wave component. The 850 mb flow (Figure 24) is nearly opposite 
to this. The upper-level GLAS GCM results also show outflow dominated by the zonal 
wind component, although the simulated meridional winds along the Equator are too 
strong. At 850 mb the GLAS GCM does simulate inflow, but contrary to the observations 
there are only very weak westerlies to the west of the reference longitude. Remarkably, the 
200 mb flow simulated by the GLA model is significantly different from both the GLAS 
model and the observations, with very weak flow along the Equator and strong meridional 
outflow regions in the subtropics of both hemispheres. The GLA model also shows a 
more realistic inflow region at 850 mb, although the meridional extent of the inflow region 
is too narrow. The UCLA GCM shows only a weak circulation more or less in the same 
sense as the observation. The model does a slightly better job for the equatorial zonal wind 
at 850 mb, but elsewhere the circulation features are not very realistic. 

The 200 mb height composite [Figure 25] for the ECMWF analyses shows a weak 
positive anomaly to the east (ahead) and a weak negative anomaly to the west (behind) of 
the oscillation center in the tropics, which are followed by enhanced subtropical response to 
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the north and south of the negative anomaly in the tropics. The wave-like response which 
appears to be emanating from the tropics is clearly identified in the extratropics The GLAS 
GCM simulates a similar wave-like pattern although the subtropical response in the 
Southern Hemisphere is not clear. The UCLA GCM also simulates a similar pattern, but 
the amplitude of the response is too weak. The extratropical response pattern in the GLA 
GCM is significantly different from that in the observations. The unrealistic subtropical 
response in this model seems to be related to the upper level divergent wind with too strong 
a meridional component in the tropics, as indicated in Figure 23. 

b. Vertical structure of the zonal component of wind 

Because of the weak signals in the UCLA model, from here on we shall focus on 
the comparison between the GLAS model and GLA model. The evolution of the observed 
vertical structure of the oscillation can be seen in longitude/height sections of the composite 
zonal wind component at the Equator for lags ranging from -10 days to +10 days (Figure 
26). The zero-lag composite shows strongest low-level inflow at about 600 mb, and 
strongest outflow at 200 mb, with the transition occurring at about 400 mb. The maximum 
low-level inflow in the u-winds is clearly above the planetary boundary layer. This 
characteristic signal is generally confined to a longitudinal strip extending from -60° (60° 
west of the wave center) to +90° (90° east of the wave center), although the upper level 
westerlies extend to +120°. The low-level westerlies penetrate the center of oscillation at 
the lowest layer. A weak westward tilt with height is apparent in the low-level inflow. 
The composites for 5 days previous to and 5 days after the arrival of the wave have a 
structure which is generally similar to but weaker than the zero-lag composite. The 
composites for +10 days and -10 days (which are nearly the opposite of each other) retain 
the upper-level structure (but shifted in longitude) of the zero-lag composite with a 
relatively weak low-level pattern. The evolution of the waves appears to have different 
stages characterized by different phase speeds. Using the upper-level easterlies as 
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reference the phase speeds are estimated to be about 10 - 12 m s 1 between ± 10 day and 
±5 day and 1 - 2 m s 1 between zero day and ±5 day. The westward tilt of zonal wind 
perturbations with height represents an eastward propagating unstable mode. It can be seen 
that the westward tilt of the wave axis becomes large as the lead or lag increases and the tilt 
is very small as the wave becomes concentrated and most intense at zero lag. This is 
consistent with the theoretical expectation of these waves as they intensify and slow down 
(Lau and Peng, 1987). Murakami and Nakazawa (1984) also observed the eastward 
propagating zonal wind perturbations with the westward tilt of wave axis. 

The GLAS GCM zero-lag composite (Figure 27) shows a realistic upper level 
configuration and a low-level easterly wind sector to the east of the wave center. The 
longitudinal extent of the pattern is somewhat less than is observed. The low-level westerly 
winds are too weak, as discussed above, although the westward tilt with height is well 
simulated in the upper troposphere above the 700 mb level. The same is true of the ± 5 
day composites; the fields are realistic except that the observed westerly inflow from the 
west is nearly absent in the model. The low-level convergence is well above the planetary 
boundary layer, as is observed. The wind pattern is quasi-stationary between -5 day and 
+5 day. There appears to be faster eastward propagation for the far-field features between 
- day and ±5 day. These are qualitatively in agreement with the observations. 

While the u-wind zero-lag composite from the GLA GCM (Figure 28) also 
indicates convergence below and divergence aloft, the convergent inflow is too strong and 
extends too high in the model atmosphere, with the maximum inflow occurring near 400 
mb and the transition from inflow to outflow above 300 mb. The 200 mb winds are also 
too weak, as we have emphasized previously. The whole pattern is slightly broader in east- 
west extent than in the GLAS model, but is not as broad as in the analyses. There is also a 
strong center of westerlies at about 350 mb, located at +120°, to the east of the easterly 
inflow branch. While westerlies at this location are in fact observed, the strong divergence 
apparent in the GLA results near +80° at 350 mb is not. The composites for lags of +5 
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days and -5 days are very similar to that at zero lag. The westward tilt of phase with height 
is very pronounced at zero lag and 5-day lag. The phase speed of the wave is quite 
uniform at about 3-4 m s 1 - 

c. Vertical structure of moisture 

The Indian Oc ean - western Pacific composites of specific humidity from the 
ECMWF analyses are shown in Figure 29. At low levels the moisture is greatest to the east 
of (ahead of) the wave center at zero lag, while the driest air is behind the wave. The moist 
air extends all the way to the surface. Because the pattern tilts westward with height, the 
moistest air lies almost at the wave center at mid-levels. The composites for various time 
lags show basically similar patterns, although translated in the zonal direction. The 
relatively slow motion of the wave in the build-up phase (between - 5 days and zero lag) 
compared to the subsequent large 5-day shift is an evidence of the wave acceleration 
eastward of the region of maximum convection mentioned in Section III. While the 
compositing technique averages the structure of the wave across the very wide Indian 
Ocean / western Pacific sector, a remnant of this dramatic acceleration in the eastern portion 
of the sector survives the averaging process. The overall signature of the oscillation is even 
broader in the moisture field (around 200° of longitude at zero lag) than was the case for 
the u-wind. The above features suggest that there is a dramatic difference in the basic 
physics between the build-up and decay phases of these intraseasonal oscillations. 

The large-scale structure of the moisture field in the GLAS GCM composite at zero 
lag (Figure 30) is generally realistic, although it is far noisier than in the observations, and 
more narrowly confined in longitude. The low-level moisture maximum ahead of the wave 
center, the drying behind it, and the westward tilt with height are all present, although the 
region of highest moisture content (over 0.4 g /kg) does not extend to the surface, as it 
does in the ECMWF analyses. While the - 5-day lag composite correctly shows little shift 
in the position of the maximum from the zero lag composite, the shift seen in the 
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observations at + 5 days is absent in the model. Instead we see a new area of high moisture 
well ahead of the wave center. Smoothing this simulated field would give a result more 
comparable to the analyses at + 5 days. 

The equatorial specific humidity composite at zero lag for the GLA GCM (Figure 
31) is similar to that from the GLAS GCM down to about 850 mb, but near the surface the 
GLA GCM shows substantial drying. Apparently, moisture convergence in the boundary 
layer itself is not required to sustain the oscillation; convergence above the boundary layer 
is sufficient. The configuration of the moisture field in the longitudinal range of -60° to 
+60° is like that of the GLAS GCM, only shifted up in the vertical. This upward 
displacement of the whole wave pattern was also seen in the u-winds discussed above. 
The rapid movement of the low-level moist tongue seen in the ECMWF composites in 
going from 0 to 5 days is also reproduced in this model. 

In general, the vertical structure of moisture appears to be closely coupled with the 
zonal circulation in the equatorial plane, in both the observations and the GCMs. The zonal 
circulation supplies moisture to the east of the oscillation center in the convergence zone 
and leaves the dry zone behind the low-level westerlies, consistent with the presence of the 
evaporation-wind feedback mechanism for maintaining the oscillation. However, it is not 
clear that this mechanism is essential for the existence of the oscillation. For example, the 
composites of the oscillation in the study of Lau et al. (1988) with and without 
evaporation-wind feedback look nearly identical. It must be stressed that determining the 
cause of the oscillation from composited fields may be quite difficult, since the relevant 
maintenance terms in the thermodynamic and moisture budgets are likely to be very small 
residuals. 
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V. TROPICAL/EXTRATROPICAL INTERACTIONS 


The connection between the tropical convection fluctuations and midlatitude 
circulation anomalies on a 30-60 day time scale has been examined in a number of studies 
(e.g. Weickmann, et al., 1985; Lau and Phillips, 1986; Knutson and Weickmann, 1987; 
Nitta, 1987, Schubert and Park, 1990). In this section we examine the global features of 
anomalous circulation in association with the tropical oscillations by one-point correlation 
using the 200 mb zonal component of wind from the ECMWF global analyses and GCMs. 
The correlations were computed between the filtered zonal component of wind and a grid 
point of the filtered 200 mb velocity potential located at the Equator and 120°E for the 
Northern Hemisphere winter (October-March) and the Northern Hemisphere summer 
(April-September), separately. The 6-month data period was used to increase the 
significance of the correlation. The 200 mb zonal wind for the GCMs and analyses are 
shown in Figures 32 and 33 for the Northern Hemisphere winter and the Northern 
Hemisphere summer, respectively. The observations show the easterly maximum in the 
Indonesian region and wide westerly zone in the eastern tropical Pacific during the 
Northern Hemisphere winter. The tropical easterlies are intensified in the Indian Ocean 
during the Northern Hemisphere summer. The easterlies in the models are somewhat too 
weak in both seasons. None of the models properly simulate the northeast tilt of the North 
Atlantic subtropical jet axis in the Northern Hemisphere winter. It is noted that the 
Southern Hemispheric subtropical jet is too strong in the GLA GCM, particularly in the 
Northern Hemisphere winter. 

Figure 34 shows the cross correlation patterns of the 200 mb zonal wind with 
respect to the 200 mb velocity potential located at the Equator and 120°E during the 
Northern Hemisphere winter from the observations. The patterns are shown separately for 
the Eastern Hemisphere in the left-hand side and for the Western Hemisphere in the right- 
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hand side for the lags of -5, 0 and +5 days. The dominant pattern is an eastward 
propagating wavenumber one structure in the tropics. The most coherent extratropical 
signal is found in the subtropical western Pacific near the East Asian jet region which 
shows a dipole pattern straddling the jet. This pattern shows an eastward shift of about 20° 
longitude in 10 days. The dipole patterns seem to be reminiscent of standing oscillations of 
global circulation modulated by eastward propagating tropical oscillations in the Indian 
Ocean and western Pacific regions. Although the extratropical signal is generally weaker in 
the Southern Hemisphere, the positive correlations in the Southeast Australia and South 
America regions appear to be related to the tropical oscillations. 

The Northern Hemisphere winter pattern of upper level wind from the observations 
[Figure 32] indicates that the western tropical Pacific is characterized by strong easterlies 
and the tropical Atlantic and eastern Pacific regions are covered by westerlies. The 
localized dipole structure of the correlation in the western Pacific may be due to the 
restricted meridional wave propagation in the easterly zone. On the other hand, the wave- 
like pattern in the Atlantic may be related to increased meridional wave propagation 
possibly in the westerly region. This is in agreement with the results of Lau et al. (1990). 

The Northern Hemisphere summer correlation patterns from the observations are 
shown in Figure 35. The tropics again show an eastward propagating wavenumber one 
pattern. The most pronounced extratropical response is found in the western Pacific and 
East Asian regions at day -5. The correlation patterns suggest that the divergent mode of 
tropical oscillations is related to anomalous anticyclonic/cyclonic circulation in the 
midlatitudes of continental East Asia. It is noted that there are some indications of 
northward propagating disturbances in the eastern Pacific. The structure of the correlation 
in the North Atlantic suggests that wave propagation from the high latitudes into the tropics 
is particularly strong 5-days prior to the extreme divergence/convergence over Indonesia. 
The connection between the Pacific dipole pattern and the pattern over the eastern Atlantic is 
not immediately clear. The wave-like patterns are not as clear as those in the Northern 
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Hemisphere winter. The extratropical signal becomes weaker in the summer hemisphere as 
the oscillation center moves into the western Pacific. The signal is pronounced near the 
Australian jet region (although not clearly seen in this projection) at day +5 which seems to 
be related to the anomalous tropical circulation south of Indonesia. We found little 
evidence for wave propagations in the Southern Hemisphere extratropics. 

The correlation patterns for the GLAS GCM are shown in Figures 36 and 37 for the 
Northern Hemisphere winter and the Northern Hemisphere summer, respectively. These 
figures show that the large-scale eastward travelling wavenumber one pattern found in the 
ECMWF analyses is not well established in the model tropics for both seasons. Since the 
realistic eastward propagating features were detected in the GLAS model from the 
irrotational part of wind in Figure 12, the stationary features of wind associated with 
eastward travelling tropical oscillations in these correlation patterns are probably associated 
with the nondivergent part of wind dominant in the model tropics. The most pronounced 
tropical/extra tropical teleconnection is found in the Australian and East Asian sector during 
the Northern Hemisphere winter. The extratropical response in East Asia propagates 
northward and appears to propagate into the Western Hemisphere. The basic state 
circulations [Figure 32 1 with westerly bias in the model tropics may provide unrealistically 
favorable conditions for wave propagation in the western Pacific. It is also indicated that 
the patterns are somewhat shifted northward compared to the observations shown in Figure 
34. The Northern Hemisphere winter patterns of upper level wind [Figure 32] indicate that 
the East Asian jet is shifted farther north in the GLAS GCM compared to the observations. 
The model simulates realistic extratropical signals in the Southeast Australia and South 
America regions. In Figure 37, the tropics and extratropics relationship is not clear in the 
model during the Northern Hemisphere summer. The model results of Simmons et al. 
(1983) and Branstator (1983) suggested that the basic state circulation is an important factor 
in determining the extratropical circulation anomaly responding to the tropical heating 
fluctuations. The substantially different basic state circulations [Figure 33] in the model 
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and the observations suggest one reason for the unrealistic response of extratropical 
circulations to the tropical oscillations. 

The correlation patterns for the UCLA GCM [Figures 38 and 39] are unrealistic, 
showing little coherence in the large scale. As described earlier, the model has much 
weaker-than-observed power in the period range of the 30-60 days. This deficiency may 
be responsible for the overall failure of the proper tropical response and the extratropical 
teleconnection in the UCLA GCM. The correlation patterns for the GLA GCM show more 
coherence but show little connection with the patterns from the ECMWF analyses. The 
tropics/extratropics interactions in the GLA GCM [Figures 40 and 41] are also not realistic. 
The unrealistic correlation patterns may be due to the inappropriate upper level divergent 
flow associated with the tropical oscillations. As indicated by the composite wind patterns 
[Figure 23], the upper level divergent flow is dominated by meridional component of wind 
in the tropics. 
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VI. DISCUSSION AND CONCLUSIONS 

A summary of .he attributes of the simulation of tropical intraseasonal oscillations in 
the 30 to 60 day range by the various GCMs is given in Table 1. Several features of 
interest emerge. Neither strong variability in the 30 to 60 day band nor an accurate 
geographical placement of such variability is enough to guarantee realistic propagation 
characteristics, as witnessed by die comparison of the GLAS and GLA GCMs. The former 
model does a reasonable job of simulating the propagation characteristics of the eastward 
propagating waves as they pass the convecively active regions in the Pacific, in spite of 
having a rather weak signal. The latter model, although exhibiting a stronger signal, tends 
to propagate it more towards higher latitudes after the wave passes the convective area. 
This appears to be related to the excessive meridional flow a. upper levels for the GLA 
model. The placemen, of the convective areas (synonymous with the region of larges, 
variability) is tied to the region where the GCM tends to produce maximum vertically 

integrated diabatic heating (precipitation) in general. 

The overall vertical structure of the oscillation, as seen in the moving composites, is 

basically realistic in both the GLAS and GLA GCMs. The GLAS model does a slightly 
better job in capturing the different stages of development of the oscillation. The dry 
conditions near the center of oscillation in the bottom layer of the GLAS GCM and GLA 
GCM are possibly related to insufficient moisture supply near the surface due to excessive 
drying of the lower troposphere by deep convection and/or a deficient wind-evaporation 
coupling at the surface. Particularly, the GLA GCM compensates for a lack of convergence 
in the planetary boundary layer by displacing the moisture convergence upward which may 
be related to the relatively weak aonal flow and an excessive meridional flow at 200 mb. 
This suggests the need to improve the boundary layer of these models, with possibly the 

inclusion of shallow convection. 
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There is a close relationship between the dominant periods that emerge from the 
spectrum and the overall vertical structure of diabatic heating associated with the 
intraseasonal oscillations. This relationship, which associates deeper heating with more 
rapidly moving waves, can be understood from earlier theory and simulations. It explains 
why the UCLA GCM, with its very deep, nearly uniform heating profile shows oscillations 
in the 15 day range instead of the observed 30-60 day range, and also helps to understand 
how the subtle difference between the GLAS and GLA GCM heating profiles (the former 
extends further downwards from its 450 mb maximum) leads to a more realistic (slower) 
phase speed for the disturbances in the GLAS GCM. This difference in heating profile is 
fairly robust, also appearing in composites of the diabatic heating in the moving frame of 
reference (not shown). The filtering which retains only 30-60 day periods is partially 
responsible for the different structures of low-frequency oscillation in the GLAS GCM and 
GLA GCM. It should be pointed out that the unrealistic heating profile in the UCLA GCM 
is not necessarily a result of deficiencies in model physics per se, and that the mid- 
tropospheric heating maxima in the GLAS GCM and GLA GCM were produced by a 
constraint inhibiting very deep convection. 

As we have emphasized, the GLAS and GLA GCMs differ from each other only in 
the numerical methods employed in solving the equations; the formulation of the physical 
parametenzations is identical in the two models. It would require further experimentation 
to determine unambiguously the mechanisms by which the differences in numerical 
techniques result in the observed differences of the tropical intraseasonal oscillations. 
Although we are not able to perform these experiments, our results strongly suggest the 
diabatic heating as one important mechanism. In particular, the geographical location of the 
dominant tropical cumulus convection and its vertical profile are important in defining the 
regions of enhanced intraseasonal variability and its vertical structure. While the tropical 
mean flows of these two models do exhibit differences, the consequences of these 
differences to the intraseasonal oscillations we present here are likely to be secondaiy to the 
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effects of diabatic heating. However, the extent of the easterlies in the time-mean flow may 
significandy affect the nature and strength of the tropical-extratropical teleconnections (e.g. 
Schubert and Park, 1990). The different basic state circulations in the model and 
observations appear to be responsible for the the unrealistic extratropical circulation 
anomalies in the GCMs. Particularly, in the GLA GCM the overall failure of simulating the 
realistic tropical-extratropical interactions appears to be related to the inappropriate upper 
level divergent flow with too strong a meridional component in the tropics associated with 
the oscillation. 

The models examined here have strongly preferred zones for diabatic heating, and 
the oscillation is forced as the region of large-scale ascent enters these preferred convective 
regions. Such a turn-on heating is observed in the intraseasonal oscillations produced by 
other GCMs with realistic lower boundaries. Finally, the present analysis is focused on the 
convective regime of ihe 30-60 day oscillations. Observations indicate that there is also a 
dry regime where the circulations seem to be maintained (with faster phase speed) in the 
absence of convection. The reasons for the fixed nature of the heating in realistic models, 
as opposed to the highly mobile convection in models with idealized lower boundaries, and 
the maintenance of the circulation in the absence of convection during part of its journey 
around the globe are major concerns for future investigation. 
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Table 1. Summary of GCM Performance 


GLAS GCM UCLA GCM GLA GCM 


Spectral Signal 

Strength weak 

Dominant periods realistic 

(30-60 days) 

Spatial Signal 

Pacific Ocean shifted eastward 

moderately weak 

Indian Ocean shifted westward 

very weak 

Local Propagation 
(Indonesian region) 

realistic 


850 mb flow 
(along the Equator) 

westerly flow west 
west of reference 
longitude weak 

200 mb flow 
(along the Equator) 

reasonable 


Vertical profile 

of heating center above 500 mb 
extends to low levels 


very weak 

moderately weak 

far too short 
(15 days) 

slightly too short 
(20-60 days) 

? 

very weak 

realistic location 
strong 

? 

very weak 

realistic location 
moderately weak 

too stationary 

disperses to high latitudes 
10 days after passing 
Indonesia 

reasonable 
but noisy 

realistic 

no signal 

zonal flow much too weak 
meridional flow too strong 

very deep 

center above 500 mb 
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APPENDIX : SEASONAL MEAN PATTERNS 
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A. Longitude-latitude maps 


A1 

A2 
A3 
A4 
A5 
A6 
A7 
A8 
A9 
A10 
All 
A12 
A13 
A 14 
A15 
A16 
All 
A18 
A19 
A20 
A21 

A22 

A23 
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Seasonal mean u-wind at 200 mb for ECMWF analyses. Contour interval is 2.0 m 
s 1 . Positive regions are shaded. 

Same as Al, except for GLAS GCM. 

Same as Al, except for UCLA GCM. 

Same as Al, except for GLA GCM. 

Same as Al, except for 200 mb v-wind. 

Same as A5, except for GLAS GCM. 

Same as A5, except for UCLA GCM. 

Same as A5, except for GLA GCM. 

Same as Al, except for 850 mb u-wind. 

Same as A 9, except for GLAS GCM. 

Same as A9, except for UCLA GCM. 

Same as A9, except for GLA GCM. 

Same as Al, except for 850 mb v-wind. 

Same as A 13, except for GLAS GCM. 

Same as A 13, except for UCLA GCM. 

Same as A 13, except for GLA GCM. 

Same as Al, except for 200 mb divergence. Contour interval is 2.0xl0 6 s 1 . 

Same as A 17, except for GLAS GCM. 

Same as A 17, except for UCLA GCM. 

Same as A 17, except for GLA GCM. 

Same as Al, except for 200 mb velocity potential. Contour interval is 2.0x1 0 6 m 2 

S' 1 . 

Same as A21, except for GLAS GCM. 

Same as A21, except for UCLA GCM. 



A24 

A25 


A26 

A27 


Same as A21, except for GLA GCM. 

Seasonal mean precipitation for GLAS GCM. Contour interval is 5.0 mm day 
Values greater than 5.0 mm day' 1 are shaded and values greater than 10.0 mm day 1 

with darker shading. 

Same as A25, except for UCLA GCM. 

Same as A25, except for GLA GCM. 
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B. Latitude-pressure sections 

B 1 Latitude-pressure sections of seasonal mean u-wind for ECMWF analyses. 

Contour interval is 5.0 m s _1 . Easterlies are shaded. 

32 Same as B 1 , except for GLAS GCM. 

B3 Same as Bl, except for UCLA GCM. 

B4 S ame as B 1 , except for GLA GCM. 
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C. Longitude-pressure sections 

Cl Longitude-pressure sections of seasonal mean diabatic heating in the tropics (10°S- 

10°N) for GLAS GCM. Contour interval is 1.0 °K day 1 . Positive heatings are 
shaded and values greater than 3.0 °K day 1 with darker shading. 

C2 Same as C 1 , except for UCLA GCM. 

C3 Same as Cl , except for GLA GCM. 
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